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Homoploid Hybridization between Native Salix cavaleriei 
and Exotic Salix matsudana ( Salicaceae) 


WU Jie’*?*, WANG Dong-chao'** , YANG Yong-ping'**** , CHEN Jia-hui'*** 
(1 Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences , 
Kunming 650201, China; 2 Germplasm Bank of Wild Species in Southwest China, Kunming Institute of Botany , 
Chinese Academy of Sciences, Kunming 650201, China; 3 Institute of Tibetan Plateau Research at Kunming , 
Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, China; 
4 University of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract; Natural hybridization between Salix matsudana and Salix cavaleriei was investigated based on populations 
from Yunnan, China. Field observations revealed that the putative hybrid, S. x heteromera had intermediate morpho- 
logies between S. matsudana and S. cavaleriei. This was further confirmed by principal component analysis. Sequence 
data of nuclear rDNA internal transcribed spacer region showed both intraspecific and intragenomic polymorphisms in 
all the three species, and S. x heteromera showed a strong additive pattern between its suspected progenitors at all 
nucleotide sites of the genotypes identified. Therefore, S. x heteromera was confirmed to be a natural hybrid between 
S. cavaleriei and S. matsudana. Flow cytometry analysis indicated that all the three species are tetraploid, and the hy- 
bridization was homoploid. Sequence data from four chloroplast datasets indicated that the hybridization was asym- 
metric, with S. cavaleriei as the maternal parent. The hybridization between the exotic common species S. matsudana 


and native rare species S. cavaleriei might increase the risk of endangerment and even extinction, indicating that the 
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introduction of Salix species should be made very cautiously. 


Key words: Salix; Homoploid hybridization; Asymmetry; Molecular evidence 


Natural hybridization is prevalent and plays an 
important role in the evolution of plants; the possible 
outcomes of hybridization include breakdown of iso- 
lating barriers; introgression; increased genetic di- 
versity; and the origin of adaptations, ecotypes, and 
species (Coyne and Orr, 2004; Soltis and Soltis, 
2009; Abbott et al., 2013). Natural hybrids usually 
show mosaic morphological characters of their paren- 
tal species; therefore, the intermediate of parental 
characters in morphology can be used to identify nat- 
ural hybrids. However, this may not always be true, 
because not all morphological characters have genet- 
ic basis. Further, when there is introgression, a hy- 
brid may be similar to one of its parental species and 
therefore difficult to identify (Rieseberg et al., 1993; 
Rieseberg and Wendel, 1993). Besides, some mor- 
phological intermediates might form via convergent 
evolution (Rieseberg and Wendel, 1993; Rieseberg 
et al., 1999). Molecular studies have shown that in- 
terspecific hybridization is more prevalent than that 
indicated by morphological and cytogenetic evi- 
dence, as reviewed by Rieseberg (1997) and Ar- 
nold (1997). Many natural hybrid species have 
been confirmed by molecular investigations, and nu- 
merous historical hybridizations have also been re- 
vealed (e. g., Hardig et al., 2000; Kaplan and Fe- 
hrer, 2007; Zha et al., 2008). 

The genus Salix L., collectively known as wil- 
lows, is a well-known taxonomically difficult plant 
taxon that consists of some 460-520 species world- 
wide, which are mainly distributed in the north tem- 
perate areas. Willows have high economic value; 
species of this genus can be used in ornamentals, 
fuel, and medicines, and are good sources of energy 
biomass as well (Fang et al., 1999; Skvortsov, 1999; 
Argus, 2010). Salix is taxonomically difficult be- 
cause of common natural hybridization, simple flow- 
ers that seldom present stable reproductive traits, di- 


oecism, and large phenotypic variation ( Rechinger, 


1992; Skvortsov, 1999; Argus, 2010). As reviewed 
by Argus (2010), there are about 120 Salix hybrids 
that have been recognized in the North American flo- 
ra (113 native Salix species are recorded in North 
American flora), and about half of these are rela- 
tively common. Indigenous species hybridize not only 
with each other, but also with introduced willows 
species. For example, the introduced Old World 
species Salix alba L. is documented to form natural 
hybrids with indigenous species Salix lucida Muhlen- 
berg and Salix nigra Marshall in New World ( Ar- 
gus, 2010). Despite the prevalent natural hybridiza- 
tion in Salix, most Salix hybrids were identified by 
morphological evidence, which might not be reliable 
as mentioned above, and seldom have been con- 
firmed by molecular evidence. A morphological and 
molecular study by Hardig et al. (2000) revealed 
that about one-third plants originally identified as 
Salix eriocephala were possible introgressants. An a- 
symmetrical natural hybridization of Populus, a 
closely related genus of Salix, was identified by 
Hamzeh et al. (2007). 

China is rich in Salix species, with about 275, 
having been recorded (Fang et al., 1999); some of 
the species described in Flora of China might be nat- 
ural hybrids. In our previous study, we found that 
Salix heteromera Handel-Mazzetti, a tree willow dis- 
tributed in limited areas of Yunnan province, China, 
always and almost only coexists with other two tree 
willows, the invasive Salix matsudana Koidzumi and 
the indigenous S. cavaleriei H. Léveillé under natural 
conditions (Fig. 1). Moreover, S. x heteromera is mor- 
phologically (e. g., leaf morphology, stamen num- 
ber, ovary stipe) intermediate between S. matsudana 
and S. cavaleriei ( Table 1, Fig.2). Therefore, we 
suspected that S. x heteromera might be a natural hy- 
brid between S. matsudana and S. cavaleriei. 

In this study, we used morphological and molec- 


ular methods to elucidate whether Salix x heteromera 
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S. matsudana 


Fig. 1 Distribution of Salix x heteromera , S. cavaleriei, and 


S. matsudana in and around Yunnan province, China 


is a natural hybrid between S. matsudana and S. cava- 
leriei, the direction of the possible hybridization and 


its impact on natural populations of these species. 


1 Materials and methods 
1.1 Plant material 

We selected five samples of the putative hybrid 
Salix X heteromera and three samples each of its sus- 
pected parents S. cavaleriei and S. matsudana from 
Lashihai, Lijiang, Yunnan, China (26°53'51" N, 
100° 7' 56" E) for sequencing of nuclear rDNA 
(nrDNA ) internal transcribed spacer (ITS) and cp- 
DNA rbcL, matK, trnD-T, and atpB-rbcL regions. 


One sample each of the above species was used for 








Table 1 Morphological comparison of the putative Salix X heteromera with the suspected parents S. cavaleriei and S. matsudana 
Characters 
Taxa Flowering time 
Leaf (length x width) /cm Number of stamens Gynophore 
S. cavaleriei 4-11 x 2-4 6-8 long March to the end of June 
S. X heteromera 5-7 x 1.2-1.4 2-5 short March to April 


S. matsudana 





1. 5-3 x 0. 6-0. 8 


none March to April 


S. matsudana 


Fig. 2 Leaf morphologyof S. cavaleriei, S. X heteromera and S. matsudana 
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flow cytometry analysis. In all, 107 specimens from 
four localities were used for morphological analysis, 
i.e., Lashihai, Heilongtan (26°52'54” N, 100°14’ 
1” E), Suhe (26°47'29”" N, 99°48'58” E), all in 
Lijiang, and Xizhou (25°51'14" N, 100°8’ E) in 
Dali, Yunnan province, China (see Table 2 for de- 
tails). All voucher specimens are deposited at the 
Herbarium of Kunming Institute of Botany, Chinese 
Academy of Sciences (KUN). 

1.2 Morphological Analysis 

Principal components analysis (PCA) of leaf 
morphology was based on five leaf characters, i. e., 
maximum blade length ( L), maximum blade width 
(W), petiole length (PL), blade length from the 
base to the point of maximum width (BL), and 
blade length-width ratio (LWR). At least five ma- 
ture leaves of each specimen were measured, and 
the average values were used in PCA. 

The PCA included data standardized for each 
trait and was performed using a correlation matrix 
without rotation; factor axes that described less than 
11% of the overall variation in leaf morphology were 
excluded from the analysis. The results were ana- 
lyzed using software SPSS 11.5 (SPSS, 2002). 

1.3 Molecular analysis 
1.3.1 DNA extraction, PCR amplification, and se- 


Total DNAs were isolated using the cetyltrimeth- 
ylammonium bromide (CTAB) method of Saghai-Ma- 
roof et al. (1984) as modified by Doyle and Doyle 
(1987). The nrDNA ITS regions were amplified by 
polymerase chain reaction (PCR) using primers 
“TTS-a” and “ITS-d” (Leskinen and Alstrom-Rapa- 
port, 1999). The direction of hybridization was de- 
termine by amplifying partial sequences of the chlo- 
roplast irnD-T, atpB-rbcL intergenic region, and 
rbcL-matK gene by using the following primes; 
“tmD°°R” and “tmT°” for trnD-T (Demesure et 
al., 1995), “atpB-1” and “rbcL-1” for atpB-rbcL 
(Chiang et al., 1998), “1F” and “1024R” for 
rbcL (Lledo et al., 1998), and “3F_KIM f” and 
“1TR_KIM r” for matK (Janzen, 2009). 

PCR was performed using a PTC-100™ pro- 
erammable thermal cycler (MJ Research, Inc.) in a 
total volume of 25 uL containing 15 uL Power Taq 
PCR MasterMix ( BioTeke Corporation), 8.5 wb 
ddH,O, 1 uL of each primer, and 1.5 pL DNA 
template. The PCR conditions included an initial de- 
naturation for 3 min at 94 C , followed by 35 cycles 
of 30s at 94°C for template denaturation, 30 s at 
50 C for primer annealing, | min at 72 °C for exten- 
sion, and a final extension period of 10 min at 72 C. 
The PCR products were purified using the PCR 








quencing Products Purification Kit (Biotype Corporation) , 
Table 2 Samples used for sequencing and GenBank accession numbers 
GenBank accession numbers 
Taxon Voucher * 
ITS rbcL matK atpB-rbcL trnD-T 

C518 KF209139-146 KF209231 KF209254 KF209243 KF209265 

Salix cavaleriei C519 KF209147-155 KF209232 KF209255 KF209244 KF209266 
C1038 ** KF209128-138 KF209230 KF209253 KF209242 KF209264 
C1030 KF209156-165 KF209233 KF209256 KF209245 KF209267 
C1047 KF209166-174 KF209235 KF209257 KF209246 KF209268 

S. X heteromera C1048 ** KF209175-184 KF209236 KF209258 KF209247 KF209269 
C1056 KF209185-193 KF209237 KF209259 KF209248 KF209270 
C1058 KF209194-203 KF209238 KF209260 KF209249 KF209271 
C523 KF209222-229 KF209241 KF209263 KF209252 KF209274 

F J C1034 KF209204-212 KF209239 KF209261 KF209250 KF209272 

eP ts 

a C1042 ** KF209240 KF209262 KF209251 KF209273 

C1099 KF209213-221 — — z — 


* All specimens collected by Jiahui Chen in Lashihai, Lijiang, Yunnan, China, and deposited in Herbarium of Kunming Institute of Botany, 


Chinese Academy of Sciences (KUN); “ Samples used for flow cytometry 
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following the manufacturer’ s instructions. The puri- 
fied PCR products of nrDNA ITS were ligated to the 
pUC18 plasmid vector, and the recombinant plas- 
mids were cloned into competent Escherichia coli 
DH5a cells ( Biotype Corporation ). The bacteria 
that contained recombinant plasmid were sequenced 
directly, and at least 8 different clones of each sam- 
ple were used. Sequences were assembled by Ge- 
neious ( Drummond et al., 2011) and aligned with 
Muscle (Edgar, 2004) , followed by manual correc- 
tion using Geneious 5.4 (Drummond et al., 2011). 
Average sequence divergence (i.e., pairwise dis- 
tance) was estimated using Kimura’s (1980) two- 
parameter method in Mega 5.2 (Tamura et al., 
2011). 
1.3.2 Determination of ploidy using flow cytometry 
DNA ploidy level of the putative hybrid Salix x 
heteromera and its suspected progenitors S. cavaleriei 
and S.matsudana was determined by comparing 
their DNA contents with those of taxa of known 
ploidy level; we used S. gracilistyla Miquel as a ref- 
erence sample, which has been reported to be dip- 
loid with 2n = 2x = 38 (Rudyka, 1990). To avoid 
the risk of error due to instrument drift, we simulta- 
neously chopped the test samples and reference sam- 
ple. Leaves (50 mg silica gel dryed leaf) of refer- 
ence and test plants were placed in a plastic petri 
dish containing 1 000 uL WPB (0.2 mol + L` Tris 
HCl, 4 mmol - L MgCl,-6H,O, 2 mmol - L ED- 
TA Na,-2H,0O, 86 mmol + L` NaCl, 10 mmol - L”! 
sodium metabisulfite, 1% PVP-10.1% Triton X- 
100, pH 7.5) for 30 min. Next, they were chopped 
using a razor blade, passed each sample through a 
30-um filter, and added 150 uL of staining solution 
(500 pg: mL"! RNase A, 1. 12 mg- mL" PI) for 15 
min in dark. Each sample was run for 2-3 min on 
Partec CyFlow Space flow cytometer. The peak of the 
reference sample was adjusted to be located approxi- 
mately at channel 100, so that the relative ploidy of 
the unknown samples could be determined by com- 
paring the peak positions of reference sample and the 


test sample by using the following ratio ( Doležel et 


al., 2007) : 
Sample ploidy = Reference ploidy x 


mean position of the sample peak 





mean position of the reference peak 


2 Resutls 
2.1 Morphological analysis 

The result of PCA of leaf morphology ( Fig. 3) 
indicated that the putative hybrid Salix x heteromera 
is an intermediate between and separate from its sus- 
pected parents S. cavaleriei and S. matsudana along 
the first factor, which accounted for 84% of the vari- 
ance observed. 
2.2 Genotypes of ITS 

The complete ITS regions of Salix x heteromera , 
S. cavaleriei , and S. matsudana were sequenced; the 
length varied from 593 to 599 base pairs (bp), and 
the aligned length was 603 bp (Fig. 4). Both intra- 
specific and intra-individual polymorphism were de- 
tected in all the three species sequenced except for a 
specimen of S. cavaleriei (e518), which had only 
one ITS repeat type. In all, 11 ITS DNA variations 
(6 in ITS1 and 5 in ITS2 region; 2 are indels and 
the other 9 are point mutations) were recognized in 
the three species. The putative hybrid S. x hetero- 
mera showed nucleotide additivity of its suspected 
parents S. cavaleriei and S. matsudana at all variation 
sites of the 11 genotypes. Further, it showed the 
most average sequence divergence that equaled the 
sum average sequence divergence of its suspected 
parents (see Table 3 and Fig. 4 for details). 
2.3 Chloroplast haplotypes 

In all, 11 haplotypes (4, 1, 2, 4 for atpB- 
rbcL, matK, rbcL, trnD-T, respectively) were de- 
tected in the chloroplast regions sequenced; 10 of 
the haplotypes of the putative hybrid Salix x hetero- 
mera were identical to those of S. cavaleriei, and one 
haplotype (a deletion in the atpB-rbcl region) was 
exclusive to S. X heteromera (see Table 4 for de- 
tails). Therefore, the hybrid S. x heteromera had 
S. cavaleriei as the plastid donor parent, and the hy- 


bridization was unidirectional, i. e., asymmetrical. 
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2.4  Ploidy level 

Flow cytometry analysis of intact leaf nuclei in- 
dicated that all the three species were tetraploid 
(Fig. 5). The diploid standard sample ( Salix graci- 
listyla) nuclei produced a single peak that appeared 
at channel 100, with average coefficient of variation 
(CV) of 9.24%, 
three test samples were around 200 ( S. cavaleriei: 
X = 198.22, CV = 4.50%; X= 
208.22, CV=4. 19%; S. maisudana; X = 209. 57, 
CV =5. 22%). Thus, the three species were conclu- 


and the peak mean channel of the 


S. X heteromera: 


ded to be tetraploid with the chromosome number of 


2n=4x=76. Our result is consistent with the repor- 


ted ploidy of S. matsudana (Suda, 1958, 1963). 
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Fig.3 Plot of leaf characters according to the first and second factor 
scores derived from PCA (factor 1 described 84% and factor 2 


described an additional 11% of the overall variation ) 
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Fig.4 Schematic diagram of sequence alignments of the ITS region in Salix x heteromera, S. cavaleriei, and S. matsudana. Sequence names 


“ . ” 
are presented as “species name_voucher_clone number 


(cav=S. cavaleriei , 


het=Salix X heteromera , mat =S. matsudana ) 
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Fig.5 Estimation of nuclear DNA content by using flow cytometry. (A) Simultaneous analysis of nuclei isolated from standard diploid species 


(Salix gracilistyla) and from S. cavaleriei; ( B) Simultaneous analysis of nuclei isolated from standard diploid species ( S. gracilistyla) and 


from S. X heteromera; (C) Simultaneous analysis of nuclei isolated from standard diploid species ( S. gracilistyla) and from S. matsudana 
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Table 3 Nucleotide positions in the aligned internal transcribed spacer (ITS) sequences in Salix x heteromera, S. cavaleriei, and S. matsudana 


Sequence region, genotype number, and nucleotide positions 
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Table 4 Chloroplast haplotype 


Sequence region and nucleotide positions 





trnD-T 


rbcL 


matK 


atpB-rbcL 


Number of 
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3 Discussion 

The biparental inherited nrDNA ITS variation of 
the putative hybrid Salix x heteromera and its suspec- 
ted parents S. cavaleriei and S. matsudana were inves- 
tigated. Our results showed that concerted evolution is 
not completed in these three species. Both intraspecif- 
ic and intra-individual polymorphisms were detected 
for all the three species. Multiple nrDNA repeats were 
common at interspecific, intraspecific, and intraindi- 
vidual levels, arising both from organismal processes 
such as hybridization and polyploidization and by ge- 
nomic processes such as gene and chromosome seg- 
ment duplication and various forms of homologous and 
non-homologous recombination ( Alvarez and Wendel, 
2003). The putative hybrid S. x heteromera showed 
both intraspecific and intraindividual polymorphism of 
nrDNA ITS region, and all the variant nucleotide 
sites were perfectly additive (i.e., chimera) of its 
suspected parents S. matsudana and S. cavaleriei ( Ta- 
ble 3, Fig.4). Stochastic genomic processes men- 
tioned above are not likely to produce such an addi- 
tive nucleotide pattern. Moreover, the ITS sequence 
diversity of S. x heteromera (0.007) was considerably 
higher than those of S. matsudana (0.004) and 
S. cavaleriei (0. 003). Further, S. x heteromera occurs 
only where both S. cavaleriei and S. matsudana are 
present; it is a morphological intermediate between 
S. cavaleriei and S. matsudana as indicated by the 
features of leaf, stamen, ovary stipe, and PCA anal- 
ysis of morphological traits. Taken together, these 
findings suggest that S. X heteromera is a natural hy- 
brid between S. cavaleriei and S. matsudana. 

Salix matsudana and S. cavaleriei also showed 
intraspecific and intraindividual polymorphism in the 
ITS region. Considering that both species are tetra- 
ploid, as shown by our flow cytometry analysis, it is 
possible that both species are of allopolyploid ( hy- 
brid) origin and have merged and maintained both 
ITS repeat types of their progenitors. Divergent re- 
peats of ITS have been reported to be clearly main- 
tained over tens of millions of years ( Baumel et al. , 


2001; Alvarez and Wendel, 2003). 
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Flow cytometry analysis revealed that Salix x 
heteromera is tetraploid and thus is a homoploid hy- 
brid. In nature, homoploid hybrid speciation might 
be a rare phenomenon; in that, parent species must 
be closely related for the homoploid hybrid to be via- 
ble, or the differences in chromosome arrangement 
might affect mitosis ( Rieseberg, 1997; Coyne and 
Orr, 2004; Abbott and Rieseberg, 2012). The pa- 
rental species S. cavaleriei and S. matsudana belong 
to Salix subgenus Salix, and sections Wilsonia K. S. 
Hao ex C. F. Fang & A. K. Skvortsov and Salix, re- 
spectively. A molecular phylogenetic study showed 
that these two sections belong to different clades 
( Azuma et al., 2000; Chen et al., 2010), indica- 
ting that S. cavaleriei and S.matsudana are not 
closely related species. This is consistent with their 
morphological differences. For example, they showed 
differences in stamen number, which is an important 
character in systematics of Salix ( Ding, 1995a; 
Skvortsov, 1999; Argus, 2010) : S. cavaleriei has 6 
-12 stamens, and S. matsudana has 2 stamens. 
They also differed ecologically; S. matsudana is 
somewhat adapted to arid and semiarid environment 
(Skvortsov, 1999) and grows along rivers or depres- 
sions amidst sand in basin and plain areas, whereas 
S. cavaleriei is a typical riparian species, requiring 
high moisture, and only occurs around streams and 
can even grow in water. Our field observations re- 
vealed that many individuals of S. X heteromera are 
big trees and apparently viable, and the hybrids 
share their habitat closely with S. matsudana; how- 
ever, the coexistence of the hybrid and its parents 
can also be observed at one site ( e.g., along a 
stream). Moreover, the flowering time of S. x hetero- 
mera overlaps with that of its parents. These findings 
indicate that S. X heteromera is not significantly di- 
vergent in ecology from its progenitors. All ho- 
moploid hybrid species that have been documented 
thus far are ecologically divergent from their parental 
species ( Abbott and Rieseberg, 2012). In addition, 
the restricted distribution of S. x heteromera ( only 


occurs when both its parents are present) suggests 


that it is sterile or its progeny is sterile and/or invia- 
ble and therefore lacks the ability to expand beyond 
its distribution range (i. e., it is not a true species). 
However, this needs to be further investigated. 

Chloroplast DNA is usually maternally transmit- 
ted in angiosperms ( Mogensen, 1996), and se- 
quencing can be used to determine hybrid origin 
(Zhou et al., 2008). Our results from the four chlo- 
roplast sequence datasets indicated that the hybrid- 
ization is unidirectional, i.e., asymmetric, with 
Salix cavaleriei as the maternal parent. Hybridization 
tends to be unidirectional at sites where one of the 
parental species is rare, because the pollens deliv- 
ered to the rare species would consist mainly of pol- 
len from the common species ( Rieseberg, 1995; 
Zhou et al., 2008). In the habitat of our putative 
hybrid, S. cavaleriei is rare and S.matsudana is 
more abundant, and S. X heteromera shared its habi- 
tat closely with S. matsudana. Under such circum- 
stance, the rare species is usually the maternal par- 
ent of the hybrid ( Rieseberg, 1995); this would 
have been the possible reason for the asymmetric hy- 
bridization observed. 

Hybridization is also often associated with habi- 
tats that have been altered by anthropogenic disturb- 
ance ( Abbott and Rieseberg , 2012). This might be 
the case in our current study; Salix matsudana is 
not documented in Floras as native species in Yun- 
nan province (Ding, 1995b; Fang et al., 1999). 
This species has long been used as an ornamental 
plant and cultivated almost all over the temperate 
zone in the world. It is widely cultivated around 
farmlands, villages, and deforestation areas. The 
cultivated plants might escape and naturalize; in- 
deed, natural population of S. matsudana is at pres- 
ent quite common in north-central Yunnan. There- 
fore, crossing occurred between the previously allo- 
patric, common and widespread S. matsudana and 
the rare and narrowly distributed S. cavaleriei. Hy- 
bridization between common and rare species might 
have severe consequences for the rare species; if fit- 


ness of the hybrid is lower relative to either parental 
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species or even sterile (i. e., outbreeding depres- 
sion), the growth rate of the rare species may de- 
cline below that required for replacement (i. e., de- 
mographic swamping ). If, however, the hybrid is 
fertile or fitness decline is negligible, the hybrid 
tends to backcross more frequently with the common 
species and may displace the rare species (i. e., ge- 
netic assimilation) ( Rieseberg and Wendel, 1993; 
Rhymer and Simberloff, 1996; Ellstrand et al., 
1999; Wolf et al., 2001). In our case, S. cavaleriei 
did not seem to be seriously threatened by its hybrid- 
ization with the exotic S. matsudana at present, be- 
cause although S. cavaleriei is rare and highly re- 
quires a moist habitat, some of its distribution range 
is not invaded by S. matsudana (e. g., Tengchong of 
Yunnan province and south Sichuan province ). 
However, if S. matsudana continues to invade the 
distribution range of S. cavaleriei by means of human 
intervention and becomes numerically superior com- 
pared to S. cavaleriei, S. cavaleriei might become in- 
creasingly endangered or even extinct through genet- 
ic assimilation and/or outbreeding depression, re- 
gardless of the fitness of hybrids. Therefore, our 
study indicated that willows should be introduced for 
purposes such as ornamentation and afforestation 
with caution, since they may cross with indigenous 
willow species and increase the risk of rare species 


becoming extinct. 
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